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A Particle and Energy Balance Model of the 
Orificed Hollow Cathode 


Matthew T. Domonkos 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 


Abstract 

A particle and energy balance model of orificed hollow cathodes was developed to assist in cathode design. The 
model presented here is an ensemble of original work by the author and previous work by others. The processes in 
the orifice region are considered to be one of the primary drivers in determining cathode performance, since the 
current density was greatest in this volume (up to 1.6 x 10 8 A/m 2 ). The orifice model contains comparatively few 
free parameters, and its results are used to bound the free parameters for the insert model. Next, the insert region 
model is presented. The sensitivity of the results to the free parameters is assessed, and variation of the free 
parameters in the orifice dominates the calculated power consumption and plasma properties. The model 
predictions are compared to data from a low-current orificed hollow cathode. The predicted power consumption 
exceeds the experimental results. Estimates of the plasma properties in the insert region overlap Langmuir probe 
data, and the predicted orifice plasma suggests the presence of one or more double layers. Finally, the model is used 
to examine the operation of higher current cathodes. 


Nomenclature 


A efr Effective emission area on the insert (nr) 

A 0 Orifice area, nr 

A r Richardson coefficient (60 A/cm 2 -K 2 ) 

dy Inner diameter of the cathode insert (mm) 
d 0 Orifice diameter (m unless otherwise noted) 

e Electron charge (1.6 x 10“ 19 C) 

E Electron energy (eV) 

E ds Electric field in the double sheath adjacent to 
the insert (V/m) 

E x Axial electric field (V/m) 

f Escape fraction of emitted electrons 

f ex Fraction of excitations causing ionization 

Id Cathode discharge current (A) 

I e Electron current to the Langmuir probe (A) 

kys Electron current from the insert region to the 

orifice region (A) 

Ie.ori Electron current emitted from the orifice (A) 

Ij Bohm ion current at sheath boundary (A) 

I i emit Emitted ion current from the orifice (A) 

Ii, or i Orifice ion current (A) 

j Current density (A/m 2 ) 

je.back Plasma electron backstreaming flux (#/m 2 ) 

je,thm Thermionic electron current density (A/m 2 ) 

j. Ion flux, (#/m 2 ) 

J, Bohm ion current density at sheath boundary 

(A/m 2 ) 

k Boltzmann constant ( 1 .38 x 1 (T 23 J/K) 

1 0 Orifice length (mm) 


InA Coulomb logarithm 

L eff Effective emission length from the insert (m) 

m^ Electron mass (9.1 1 x 10“ 31 kg) 

m ; Ion mass (2. 1 9 x 1 0“ 25 kg for xenon) 
m Mass flow-rate (mg/s or seem) (1 seem = 
0.097 mg/s) 

M e Atomic mass of an electron (1/1836) 

M; Atomic mass of an ion (1 3 1 for xenon) 

M 0 Atomic mass of a neutral (131 for xenon) 
n x Number density of species x(cm“ 3 or nT 3 ) 

N eo Number density-rate of electron impact 

ionization events (#/m 3 -s) 

Qconv Convective power loss (W) 

fie, back Convective power loss due to backstreaming 

electrons (W) 

q e con Convective power loss due to electron current 
(W) 

fie.thm Convective power input by thermionic 

electrons (W) 

q ex Power loss due to radiation of excited states 
(W) 

fii.ioss Convective power loss due to ion loss (W) 

qi.on Convective power input by orifice ions (W) 

q i0 n Ionization power loss (W) 

qohime ohmic heat generation (W) 

r G Orifice radius (m) 

R Resistance (Ohms) 

Rs P Specific gas constant (63. 1 J/kg-K for xenon) 
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U ex Average excitation energy (eV) 

v e , v e Electron thermal speed (m/s) 

v t Average ion speed (m/s) 

V ds Voltage across a planar double sheath (V) 

V p Plasma potential (V) 

e 0 Permittivity of free space (8.85 x 1(T 12 F/m) 

(|) eff Effective work function (eV) 

([); Ionization potential (12. 12 eV for xenon) 

<|) 0 Material work function (eV) 

y Ratio of specific heats (5/3 for xenon) 

r) Plasma resistivity (Q-m) 

K Wavelength of pyrometer optics (m) 

X,d Debye length (cm or m) 

X ee Electron mean free path for self-collision (m) 
\ e x Electron mean free path for excitation (m) 

^,i on Electron mean free path for ionization (m) 

p. Permeability (H/m) 

Ho Permeability of free space (4k x 1(T 7 H/m) 
v a _ b Collision frequency of a with b (Hz) 
p Mass density (kg/m 3 ) 

p 0 Mass density in the orifice (kg/m 3 ) 

o Self collision cross section for ions and 
neutrals (nr) 

o a>b Collision cross section for a colliding with b 
(m 2 ) 

o ex Electron-impact excitation cross-section (nr) 

Ci on Electron-impact ionization cross-section (nr) 

Introduction 

The current carrying capability and long-life of orificed 
hollow cathodes has made them the primary choice for 
electric propulsion systems from 100-W to several tens 
of kilowatts. 1 For low current devices, the power and 
propellant consumption of the cathode, in addition to 
life, are the design drivers. At high current, the ion flux 
may become sufficient to limit cathode life by erosion, 
particularly for the ion thruster discharge cathode. 2 
Despite the volume of experience with hollow cathodes 
in the electric propulsion community, thorough 
understanding of the plasma phenomena within and 
near hollow cathodes remains elusive, due in part to the 
small scale of the devices. Cathode life prediction is 
hindered by the fidelity of the measurements and 
calculations of plasma properties within and near the 
cathode. Several models have been developed 
previously to examine the insert region or orifice region 
plasmas separately. A hollow cathode model which 
considers both the insert and orifice plasma regions is 
presented in this paper. The model is used to examine 
plasma properties within the hollow cathode and 
provides a starting point for examining the plasma 
downstream of the cathode. 


Model Description 

The hollow cathode model presented here is based on a 
control volume approach to the orifice and insert 
plasmas respectively. An overview of the cathode 
model is presented in Figure 1. Conservation of ions 
and energy, as well as current continuity are applied to 
each control volume. These sets of equations provide a 
basis for determining the electron temperature and 
number density in each region and the insert 
temperature in the insert region. The orifice model, 
which is based on the one originally developed by 
Mandell and Katz, 3 is presented, followed by the model 
of the insert region. 


Orifice Model 

The orifice model assumes that all of the cathode 
current is conducted through the orifice, and 
consequently ion collection to external cathode surfaces 
is neglected. The plasma within the orifice is treated as 
quasineutral, and the ions lost to the cathode and the 
insert region are replaced through electron impact 
ionization. Energy is deposited into the plasma via 
ohmic heating, and ionization and convection are the 
energy loss mechanisms considered. Mandell and 
Katz 3 developed a simple model of the orifice processes 
based on these conservation equations, and the present 
derivation includes current continuity. 


Conservation of Ions 

The insert region is modeled as a cylindrical control 
volume bounded by a cathode sheath at the orifice 
radius, by a double sheath at the entrance to the orifice, 
and by the exit plane of the orifice. Figure 1 illustrates 
the orifice model approximations. The chamfer is 
neglected in the analysis of the control volume. In the 
formulation of the model, the plasma properties are 
assumed to be uniform within the orifice, and the radial 
electric field is neglected. The axial electric field, 
however, is implicitly included in the orifice model. 


The steady-state continuity equation applied to the ions 
states that the creation and in-flux of ions is balanced 
by the out-flux. 


dn j q 

( dn. 

( dn, 'l 

( dn, \ 


— L + — H 


dt 

\ (*t ) i 0n i zat i 0n 

V dt ) ^ _ jjnx 

V dt Jout-fita 


( 1 ) 


where n; is the ion number density. The second term is 
negligible. The double-layer at the entrance to the 
orifice region repels ions, and the flux of ions from 
downstream of the orifice, calculated from measured 
plasma properties, is much smaller than either of the 
other two terms. 4 ' 3 Consequently, the term describing 
ion flux into the orifice is omitted from the model. In 
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order to evaluate the first term in Equation (1), a 
collision analysis was performed, and the derivation 
appears in Domonkos. 6 For simplicity, electron-impact 
ionization was assumed to be the only mechanism for 
creating ions within the orifice. Equation (2) describes 
the ionization rate 


f ^ n i \ 

L*Li ) 

f m ‘ )' 

[nr J . . „ 

v / ionization 

{ 4 Z «J 

(imtJ 


2 

co mv e 

x4m e n 0 ^dv e v i e a i on{v e )e 2kr ‘ 

0 

where d 0 is the orifice diameter, 1 0 is the straight 
channel length of the orifice, m,, is the mass of an 
electron, T e is the electron temperature, n e is the 
electron number density, n 0 is the neutral number 
density, v e is the electron velocity and G ion (v e ) is the 
cross section for electron-impact ionization. Equation 
(2) neglects the drift induced by double layer because in 
most cases the potential difference is small. The neutral 
number density is calculated by assuming choked flow 
at the exit of the orifice 


ions out of the control volume is usually less than ten 
percent of the total ion loss term. Consequently the 
density gradient at the exit of the orifice determines the 
ion motion. The model by Capacci, et al. 8 included a 
sheath at the downstream end of the orifice which is 
neglected in the present derivation. The sudden 
expansion of the plasma at the orifice exit approximates 
the conditions of a decaying plasma, and ion transport 
was hypothesized to occur by ambipolar diffusion. The 
transport of ions downstream of the orifice in a 
quasineutral plasma is essential to the theory of spot 
and plume-mode operation proposed by Mandell and 
Katz/’ In the present model, the emitted ion current is 
also used to estimate the emitted electron current 
beyond the keeper. 

Current Continuity 

Current continuity is described at the orifice entrance 
and at the orifice exit. At the entrance, the total 
discharge current is calculated as the sum of the 
electron current from the insert region and the ion 
current to the walls of the orifice less the ion current 
emitted into the cathode-to-keeper gap. 


n„= 7— (3) 

m 0 a 0 ^rJ 0 

where M 0 is the atomic mass of the propellant, A 0 is the 
orifice cross sectional area, and the neutral temperature 
is assumed to be equal to the orifice plate temperature. 
The orifice plate temperature was estimated based on 
the experimental data. The integral in Equation (2) is 
evaluated numerically using the experimental values for 
the cross-section by Rapp and Englander-Golden. 7 The 
ion flux out of the control volume is determined using 
the Bohm condition at sheath boundaries and the 
thermal flux at the downstream boundary of the control 
volume. 



= — = — + r ‘ )) 


= 0.61 nc i {ixrj,, +rcr‘ ) + -»v ( icr \ ) 


(4) 


« , 

= 0.61/! ( 2 !tr l + nr 

- 1 ... V 



This equation implicitly includes the effects of the 
radial electric field with the use of the Bohm criterion. 
At the exit of the orifice, the density is assumed to be 
sufficient to shield the bulk of the ions from drifting 
toward cathode surfaces, and the thermal flux of ions is 
used. The Debye lengths in the orifice indicate that the 
bulk of the orifice volume is free of the strong radial 
electric field near the orifice wall. The thermal flux of 


I D I e,ins ^ i,ori ^ i .emit 


(5) 


The orifice ion current is calculated as part of Equation 
(4), as is the emitted ion current. The insert electron 
current is determined using Equation (5). At the 
downstream end of the orifice, the discharge current is 
assumed to be the difference between the emitted 
electron and ion currents. 


I D I e.ori I i. emit (^) 

Equation (7) combines Equations (5) and (6). 

1 d = let* + Q-6\n e c t Jt{lr 0 l 0 + r 2 )~\n e v t [nr 2 ) 

j (7) 

= I e,on~-n e v i {^r 2 ) 

This type of accounting is necessitated by estimation of 
the ion current in the orifice being on the order of ten 
percent of the total current. The electron current is 
assumed to increase linearly between the orifice 
boundaries. While the rest of the orifice plasma model 
treats all of the properties as spatially invariant, the 
axial variation of the electron current influences the 
ohmic heating within the orifice. 


Conservation of Energy 

Conservation of energy within the orifice control 
volume is approximated by equating the ohmic heating 
with the losses from ionization, radiative decay of 
excited states, and convection. 
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q Ohmic dion Qex ~^Qconv (^) 

The ohmic heating is calculated using the resistivity 
based on electron-neutral and electron-ion collisions, as 
shown in Equation (9) 9 


where v is the sum of electron-ion and electron-neutral 
elastic collision frequencies. As an approximation, the 
electron-ion collision frequency is calculated from the 
formula in the NRL Plasma Formulary. 10 

_2 

v . = 3.91nA fiT 2 (10) 

ei e e ' ' 


uniform in the orifice control volume for simplicity. 
The formulation of the ohmic heating term 
accommodates the axial variation of the electron 
current. Thus for x < l c 

Ie(*) = Ie^+bx ( 14 ) 

where I ejills and b are determined by the requirement for 
current continuity, and are defined in Equations (15) 
and ( 16), respectively. 

I e , ins =I D -Q.6\n e c t {27Cr 0 l 0 +Kr;) (15) 

b = 1 - 1/ ’" ,s = y-[o.6 1 n e c, (2 jrr 0 l 0 + nr] )] (16) 


where the dimensions are consistent with those listed in 
the Nomenclature and the Coulomb logarithm is 


In A = 23-ln 




( 11 ) 


Calculation of the ohmic heating of the electrons is 
modified to account for the variable electron current 
along the length of the orifice. 


. = t 2 r = JL 

“ ohmic 1 /V j 

A„ 


Ie,Jo+Ie,^l 2 0 +-b 2 ll 


( 17 ) 


An effective value of the electron-neutral elastic 
collision cross-section is estimated by numerically 
integrating the velocity dependent elastic collision 
cross-section data from Brode over the Maxwellian 
electron population in the same manner as 
Equation (2)." The effective cross-section is then used 
to evaluate the electron-neutral elastic collision 
frequency as defined in Equation (12), 


The ionization power loss due to the flow of ions out of 
the control volume is the product of the ionization rate, 
the charge, and the ionization potential: 

q ,o„ = { ^7] e<j) t = N eo (/rr 2 l 0 ) eQ t (18) 

V / ionization 

where N eo is determined from Equation (2), and <]); is 
equal to 12.12 eV for xenon. 


K, o ={n o -n e ){(T e , 0 v e ) (12) 

where the neutral number density based on flow, no, is 
reduced by the electron density under the quasineutral 
assumption, and the velocity of the electrons is defined 
as 



The factor of V2 typically used to calculate the 
collision frequency is omitted from Equation (12) 
because the electron velocity is assumed to be much 
greater than the neutral velocity, thus defining the 
relative velocity. Since the properties within the orifice 
are assumed to be constant everywhere, Equation (9) 
defines the plasma resistivity in the orifice. The 
electron current in the orifice is assumed to increase 
linearly toward the orifice exit due to the creation of 
ions and electrons through collisions with neutrals 
along the length of the orifice. The variation in the 
electron current implies an axial variation in electron 
density, however the electron density was treated as 


The second term on the right-hand side of Equation (8) 
describes the power lost from the free electrons as they 
excite bound electrons within neutral xenon. The 
energy lost by free electrons in excitation events is 
calculated similarly to the ionization energy loss. The 
excitation collision rate in the orifice is (19) 



x4mijn 0 - n e )j dv^o^fvje 


0 

where the collision cross-section data were taken from 
Hayashi. 12 The total energy loss of the free electrons to 
the excitation of the bound electrons is 


dex = 


dn a 

dt 


eU„ 


( 20 ) 


where the average energy lost in an excitation event, 
XJ ex , is a constrained parameter. Similarly to Mandell 
and Katz, 3 the average excitation energy is set to 10 eV. 
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While this value is representative, provided the 
overwhelming majority of the excitations are from the 
ground state, the average excitation energy decreases as 
the fraction of secondary excitations increases. In order 
to rigorously evaluate the power loss due to neutral 
excitation, a detailed ionization model is required. 
Development of such an ionization model was beyond 
the scope of the present investigation, and reduced 
accuracy of the results can be tolerated when using the 
model to assess the relative capabilities of a given 
cathode design. 


Current Continuity 

In the insert region, the current continuity statement is 
written by equating the total operating current, I D , with 
the current emitted from and collected at the cathode 
surface. This approach necessitates inclusion of the ion 
current in the orifice, and the statement is therefore 
dependent upon the results of the orifice model 
calculations. The mathematical expression representing 
current continuity is written in Equation (22). 

^ D 1 e.emit ^ i.coll ^ r, coll ^ i .emit ( 22 ) 


The convective power loss considers the temperatures 
of the electrons entering the control volume from the 
insert region and those exiting the orifice: 


The first term in Equation (22) refers to the thermionic 
electron emission current from the insert which is 
determined using the Richardson-Dushman equation. 


tfconv e,ori e,ori e,ins e,ins ( 21 ) 

The electron temperature in the insert region is solved 
for in the insert model. In general, the total convection 
term is a small fraction of the energy lost in neutral 
excitation. 

The ion continuity (1) and energy equations (8) are 
solved iteratively by using a numerical goal seeking 
program which varies the electron temperature and 
number density, respectively. First the electron 
temperature is varied to satisfy Equation ( 1 ). Next, the 
number density is determined, satisfying Equation (8). 
If the solution of Equation (1) holds within established 
limits, the solution is considered to have converged. 

After the orifice model solution converges, several 
Knudsen numbers are calculated to verify that the 
plasma is indeed collisional, and the Debye ratio and 
particle count in a Debye sphere are evaluated to 
determine if the ionized gas meets the conditions for a 
plasma. The results showed that the plasma was 
collisional with Knudsen numbers generally much less 
than one for plasma-related collisions. The Knudsen 
numbers for neutral and ion self-collisions were as high 
as a few tenths. This is considered borderline for any 
continuum theory describing their properties. 

Insert Model 

In the same spirit as the orifice model, the plasma in the 
insert region is approximated by a control volume with 
uniform plasma properties throughout. Figure 1 
illustrates the basis for the insert model. The efficiency 
of the ionization processes within the insert region is 
expected to contribute to the overall performance. In 
addition to the statements of ion continuity and energy 
conservation, the insert model has an explicit statement 
of current continuity. These three equations still leave 
four free parameters, and the choice of those parameters 
is discussed. 


I e , mU =I^={2^ t L eff )A R Tle kT - ( 23 ) 

The value of the Richardson coefficient was set 
constant at 60-A/cm 2 -K 2 which is consistent with 
Goodfellow. 1 "' The work function is considered as a 
bounded free parameter. Typical values for the work 
function of the insert range from 2.0 to 2.1-eV, and the 
variation of the thermionic current density as a function 
of temperature is plotted in Figure 2. 14 An effective 
work function is used to account for the Schotky effect 
which acts to reduce the apparent work function of a 
material in the presence of a strong electric field. 



A double sheath analysis for an electron emitting 
surface estimates the electric field in Equation (24) to 
be 15 ’ 16 


E* = 



2 



1/2 


4 


( 25 ) 


The plasma potential in Equation (25) is a restricted 
parameter in this investigation. Experimental data 
showed that the plasma potential varied between 8 and 
1 2-V above the cathode potential depending upon axial 
position. 22 The effect of the plasma potential on the 
thermionic current density is depicted in Figure 3 for 
plasma parameters similar to those expected. While the 
emission current density is relatively insensitive to the 
plasma potential, the Schotky effect increases the 
emission current density significantly. 


The ion current collected at the cathode surfaces 
includes the ion flux in the orifice and insert regions. 
Ion backstreaming to the insert region is also calculated 
using the Bohrn criterion, and these ions are assumed to 
recombine on the cathode surface in the insert region. 
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Additionally, the plasma ions, generated in the insert 
region, also contribute to the ion current. In the context 
of the control volume approximation, ions enter the 
cathode sheath at the Bohm velocity and diffuse 
upstream at the thermal velocity. Ions diffusing 
upstream were also assumed to recombine on the 
cathode surface. 

Ion emission and electron backstreaming are both small 
terms, however they are included for completeness. 
Electron backstreaming to the insert surface is 
described using Equation (26) 

1 _ ~~ L 

J e,back = ~f n eV e eA eff e kT - ( 26 ) 

which is on the order of a few milliamperes typically. 
The effective surface area for electron collection, A eff , 
includes both the effective emission region of the insert 
and the orifice plate surface. Diffusion of electrons 
upstream of the insert control volume is neglected 
based on the assumption that the axial electric field 
overwhelms the tendency of any electrons to diffuse in 
that direction. 


Ion Conservation 

This species specific form of conservation of mass is 
essentially the same as that used in the orifice mode. 
Equation (1). The primary difference is the 
consideration of two electron populations: one a beam 
from the thermionic emission through a double sheath, 
and the other the Maxwellian plasma electrons. The 
contribution of the plasma electrons to the creation of 
ions and excited atoms is exactly that prescribed by 
Equations (2) and (19) using the insert effective 
volume. 


The electrons emitted from the insert are referred to 
here as the primary electrons. The primary electron 
speed distribution function is approximated as 


f(v e ) = 4x 




2nkT 


v]e 


( 27 ) 


where Vd is the electron drift velocity induced by the 
sheath voltage, and the temperature of the primary 
electron beam is assumed to be the insert temperature. 
Since the primary electrons are created at the insert. 


excitation and ionization collisions deplete the 
population with increasing distance from the insert. 
The mean free path for either of these collisions is 
calculated as shown in Equation (28) 


A = — 

m 7 


nAn\ 


2nkT, 


\dv e v]a{v e )e 




<v.) 


( 28 ) 


where the number density is that of the target. Figure 4 
depicts the electron impact cross-sections for 
xenon. 17,18,19 Below 20-V, primary electrons emitted 
from the insert are more likely to undergo an excitation 
collision than an ionization event. Measured plasma 
potentials within the insert region are 14-V and lower. 22 
Evaluation of Equation (28) for the electron-impact 
excitation cross-section and the electron-impact 
ionization cross-section indicates that the mean free 
path for the ionization interaction is at least 10 times 
greater than for excitation in the range of plasma 
potentials observed experimentally. 22 For simplicity, a 
planar geometry is used to approximate the primary 
electron number density at the mean free path for 
electron-impact ionization. 




(r 0 ) l- erf 



( 29 ) 


For X I0n »X ex , the contribution of electron-impact 
ionization by the primary electrons is negligible, and 
consequently this ionization mechanism was neglected. 
Nevertheless, the primary electrons are expected to 
contribute strongly to ionization. The only mechanisms 
considered for ion creation are electron-impact 
ionization by plasma electrons, and step-wise excitation 
by both the primary and plasma electrons, leading to 
ionization. The latter phenomenon is considered 
without a rigorous derivation of the contribution of this 
term. However, calculated excitation collision 
frequencies for the primary electrons in the insert 
region are on the order of or greater than the Einstein 
coefficients for neutral xenon, indicating the feasibility 
of a multi-step ionization process. 20 In the model, a 
fixed percentage of the excitation collisions are 
assumed to create ions. The resulting expression 
describing the ion production in the insert region is 
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= Wi L, 


eff) 


ionization 

f 


2nkT„ 


% 


j 


4m e n 0 dv ( ,v„a ion (v e )e 


- fMn 2 -in -2A ex f)i eff \ 


2kT p 


X 


2 nkT 


(30) 


x4m >e, pr (no - n e )J dv e v]i T ex (v e )e 


The percentage of excitation collisions which create 
ions f a is one of the free parameters in the insert 
model, and was expected to be a few percent at most. 
The effective volume for the excitation of neutrals is 
limited to two mean free paths from the insert. This 
adjustment accounts for the depletion of primary 
electrons beyond this distance from the insert. 


Ion creation by collisions of electrons with neutrals and 
ion in-flux from the orifice are balanced by ion 
convection across the control-volume boundaries. 
Plasma ions enter the cathode sheaths at the insert 
radius and at the orifice plate with a flux determined by 
the Bohm condition. Ion loss from the control volume 
is also calculated by the thermal flux to the upstream 
boundary. The various fluxes calculated in this section 
are also useful in determining the conservation of 
energy within the orifice. 


Conservation of Energy 

The conservation of energy in the insert region is 
written for the plasma as a whole. Energy enters or is 
released in the control volume by electron convection 
from the insert, ion convection from the orifice region, 
and Ohmic heating. The energy convected with the 
thermionic electrons is 


voltage drop across the double sheath at the boundary 
between the insert and orifice regions, V ds , was 
estimated by Capacci, et al. 8 to be 


9I D kT e 

15A 0 n e e 2 



(33) 


Ohmic heating is calculated in a similar manner as 
described for the orifice. The plasma resistivity is 
determined assuming uniform density and temperature. 
The electron current to the orifice calculated from 
Equation (5) is used for the Ohmic heating calculation. 
Since the current is carried mostly in the radial 
direction, an average cross-sectional area for current 
conduction is used. The resulting expression for the 
Ohmic heating is 


(34) 

The energy loss mechanisms from the insert region are 
ion loss, excitation collisions, electron convection of 
the current, and electron backstreaming to the insert. 
The ion loss term is calculated based on the ion flux to 
the boundaries of the control volume. 


i i, loss 


0.61^c j e(2^.I €i 


+ K ( r? - r 2 



1 _ 

+ ~n e v t e 



The energy lost in excitation collisions is calculated as 
the product of the rate of excited neutral production and 
the average excitation energy. 


<l„ = 


dn 


dt 


el T 


= k (». - ", )kO, ( 2 * r A f ) + » (» - ", )( v . ff .,) p , ] eU „ 

(36) 


q,n m =Ie Am \V f +^-] > (31) 

where the fall voltage, Vj , is approximated by the 

plasma potential. The energy convected with the 
orifice ions is 

9i,ort = (0.6 1 n e c t eMr 2 ^ [v ds + ^ j (32) 

where the average thermal energy of the ion flux to 
the upstream orifice boundary is 2kT t because the 
energy transport favors the high-energy particles. 21 The 


where the average excitation energy U ex is another free 
parameter. The value of U ex is somewhat lower than in 
the orifice because the energy lost by electrons in 
secondary excitation collisions, which play an 
important role in ionization in the insert region, is less 
than 4-eV, thereby reducing the average energy 
compared to the case where only ground-state 
excitation occurs. For most of the results presented, U ex 
is set equal to 5-eV for the insert region, and a 
sensitivity analysis was conducted. The energy lost by 
electron conduction of the current is 

q e ,con=—UT e (37) 

e 
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where T e is the Maxwellian electron temperature in the 
insert region. The insert electron temperature 
calculated from this model could be used as an input for 
the orifice model. Since the orifice model is relatively 
insensitive to this parameter, an average value of the 
insert temperature is used to expedite the solution. 
Finally, the electron backstreaming component of the 
energy equation is 



The solution method for the insert model is essentially 
the same as for the orifice model. Ion continuity is 
satisfied by varying the electron temperature. 
Conservation of energy is achieved by solving for the 
electron number density, and the insert temperature is 
varied to achieve current continuity. Both the orifice 
and insert models were nm with inputs matching those 
for the experimental set-up, and the results are 
discussed in the next section. 


Model Predictions 

The inputs for the orifice and insert models are the 
current, flow-rate, cathode material, and cathode 
geometry. By using the same inputs available to an 
experimental investigation, the results of the model are 
directly comparable to the experimental data. In this 
investigation, the geometry of a 3.2-mm diameter 
hollow cathode 22 was input to the orifice and insert 
models, and the flow-rate and current were varied over 
the range tested experimentally. Additionally, the 
models were also used to evaluate the effects of 
geometry and material changes for low-current 
cathodes. While the standard operating procedure for 
the models was to vary only the parameters controllable 
in an experimental situation, several free parameters 
remained. The model was evaluated to test the 
sensitivity of the various solutions to the free 
parameters. 


Sensitivity of the Orifice Results to the Free 


Parameters 


The orifice model neglects the creation of ions by step- 
wise excitation, and the calculated Knudsen numbers 
for this type of collision indicate that this phenomenon 
occurs infrequently. Nevertheless, the effects of a fixed 
percentage of the excitation collisions creating ions 
were evaluated, and the results are plotted in Figure 5. 
Unless stated otherwise, all the results presented here 
refer to the cathode AR6 of Reference 22. The number 
density shows the largest change, monotonically 
increasing with f ex . In terms of a figure of merit for 
cathode optimization, the ohmic heating and excitation 
energy loss both scaled with nearly a ten percent 


reduction for every five percent increase in f ex . Since 
the Knudsen numbers suggested that step-wise 
ionization was the exception in the orifice, f ex was set to 
zero for subsequent calculations. 

The average excitation energy is also considered a free 
parameter since the distribution of the excitation 
transitions was unknown. Excitation from the ground 
state in xenon costs between 8.44 and 10.40-eV. 23 If 
this event dominates the excitation collisions, then the 
average excitation energy is on the order of 10-eV. 
While the preceding paragraph discounted the 
contribution of step-wise excitation in the orifice, the 
effect of reducing the average excitation energy was 
evaluated to examine the potential consequences of this 
simplification, and the results are shown in Figure 6. 
The most sensitive parameters are the number density 
and the power consumption. The power consumption 
was expected to scale with the excitation energy. The 
number density increases as the excitation energy 
decreased, ensuring that the convective power terms 
account for the difference between the ohmic heating 
and the excitation loss. The previous discussion 
concerning the orifice Knudsen number for excitation 
supports the conclusion that ground state excitations 
dominate, and the average excitation energy was set to 
1 0-eV for all subsequent calculations. 


Sensitivity of the Insert Results to the Free 


Parameters 


By contrast, the excitation processes in the insert region 
are essential to plasma generation. The primary 
electrons undergo excitation collisions almost 
exclusively. Figure 7 depicts the sensitivity of insert 
region plasma parameters and power consumption to 
the value of f ex in the insert region. Electron convective 
power consumption decreases by ten percent for an 
increase of two percent to f ex , while the number density 
increases by twenty percent over the same interval. 
Although the number of excitation events possible in 
xenon argue that only a small fraction create ions, an 
ionization model is needed to rigorously determine the 
value of f ex . The fraction of excitations creating ions 
was set to 0.05 for all subsequent calculations, since the 
results agreed well with the experimental data. 


As important as the value of f ex is to the solution of the 
insert region parameters, the average excitation energy 
in the insert region influences the solution to a greater 
degree. Figure 8 depicts the sensitivity of the plasma 
properties and energy terms in the insert region to 
changes in the average excitation energy. While the 
plasma properties and energy terms are highly sensitive 
to U ex , the preponderance of the secondary excitations 
in the insert region is expected to reduce the average 
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excitation energy compared with the ground-state value 
used in the orifice calculations. Over the range from 3 
to 8-eV, the change in power consumption by 
excitations scales with U ex . The number density is 
highly sensitive to the average excitation energy in the 
insert region, changing by about ten percent per 
electron-Volt. Since many of the possible transitions 
between excited states in neutral xenon involve less 
than 1-eV, 23 the lack of detailed accounting of the 
population of the excited states makes the value chosen 
for U ex somewhat arbitrary, although it should be 
weighted for the ground-state transitions which enable 
secondary excitations. The average excitation energy 
was set to 5-eV in lieu of a more rigorous evaluation. 
This number accounts for electron collisions with 
excited states while a lower number requires that the 
plasma approach optically thick conditions. The 
experimental data indicated that the insert region 
plasma was far from optically thick conditions. 

The plasma potential is a free parameter in the insert 
model, although experimental data put bounds on the 
value chosen. Furthermore, the sharp axial gradient in 
the plasma potential made the choice of a volume- 
averaged value perilous. The plasma potential was 
varied from 7 to 12-eV to obtain the data in Figure 9. 
Again the number density is highly sensitive to this 
parameter as is the convected energy from the 
thermionic emission current. The latter finding is 
expected since the reduction in the work function by the 
Schotky effect scales with number density. While the 
insert region plasma potential near the insert was as 
high as 12-V, a value of 8-V was used for subsequent 
calculations both in consideration of the volume 
averaged value and because large values for V p 
degraded the stability of the model. 

The final free parameter considered here is the material 
work function of the insert. While the available data 
confine the work function to between 2.0 and 2.1-eV, 
this range strongly affects the requisite cathode 
operating temperature as shown in Figure 10. A change 
of 0.1 -eV in the work function leads to more than a 
50-K increase in the insert temperature. It should be 
noted that in a real system, the plasma potential 
depends upon the material work function, and the 
variation of either of these parameters independently is 
artificial. A value of 2.00-eV was chosen for 
subsequent calculations primarily due to an apparent 
over-prediction of power consumption; a low-work 
function reduces the power consumption gap at high 
flow-rates. 

Sensitivity of Predicted Cathode Power 

A comparison of the sensitivity of cathode power 
consumption to variations in several of the free 


parameters is depicted in Figure 1 1 . Cathode power is 
most strongly tied to the choice for the average energy 
of excitation collisions within the orifice. The fraction 
of excitation collisions resulting in ionization in the 
orifice and the plasma potential in the insert region 
exert comparable influence on the predicted cathode 
power consumption. The average energy transferred in 
an excitation collision in the insert volume influences 
cathode power consumption in a non-linear fashion, and 
the magnitude of dependence of cathode power on the 
average excitation energy in the insert region is less 
than the dependence upon the other free parameters. 
Plots similar to Figure 1 1 provide insight into the 
suitability of the model predictions for use in cathode 
design. If the variability in the free parameters is 
minimized, the model predictions become more 
reliable. 

Comparison of Cathode Power Consumption with 
Experimental Data 

The model predictions are directly compared with 
experimental data for a 3.2-mm diameter hollow 
cathode in Figure 12. The average excitation energy 
and the fraction of excitation collisions resulting in 
ionization were varied in an attempt to match the 
experimental data. The predicted power consumed in 
the orifice is plotted in Figure 12, along with the total 
power consumed experimentally. In order for the 
predicted power in the orifice to be consistent with the 
total power measured, the free parameters were set well 
beyond the bounds discussed above. The control 
volume approach appears to be too simplistic to capture 
the subtleties of hollow cathode operation at low 
current. Observation of Figure 12 indicates that in spot 
mode, the measured cathode power consumption was 
nearly invariant with flow rate. The positive slope of 
the predicted power consumption results from 
Equations (18) and (19), where the increase in neutral 
density from the flow rate results in additional power 
losses for the control volume. 

Insights from Calculated Plasma Parameters 

The ability of the model to predict the plasma 
conditions within the cathode was also evaluated. 
Figure 1 3 depicts the variation of the plasma properties 
within the cathode AR3 from Reference 22. The insert 
region is the only location where the experimental data 
overlaps with the model. In the insert region the 
predicted, bulk-averaged electron number density and 
temperature match the experimental results. Since 
models are currently the only means to evaluate the 
plasma properties within the orifice and the model 
appears to predict plasma properties to within a small 
error, the cathode model presented here was used to 
gain insight into the processes governing current 
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conduction in the cathode. The large gradient in 
number density and the modest gradients in electron 
temperature both upstream and downstream of the 
orifice suggest the presence of at least one double layer. 
The model assumes a double layer at the orifice 
entrance, and the rapid expansion at the exit of the 
orifice provides evidence of a double layer in this 
region as well. The downstream double layer may 
enable the electrons in the orifice to gain sufficient 
energy to ionize local to the orifice, thereby creating the 
spot and facilitating current conduction. 

Comparison with the Space Station Plasma 
Contactor 

Although the cathode model was originally developed 
as a design tool for low-current hollow cathodes, an 
attempt was made to assess its performance in 
predicting the operation of existing hollow cathodes. 
The Space Station Plasma Contactor (SSPC) is one of 
the most widely tested hollow cathodes and is ideal for 
comparison purposes. Figure 14 shows the predicted 
orifice power consumption in the SSPC, and the total 
power consumption of the SSPC at 3.0-A discharge is 
also shown for comparison. 24 Unlike the low-current 
hollow cathode results, the predicted orifice power is 
less than the total experimental value. Additionally the 
predicted trends of orifice power consumption with 
flow-rate more closely approximate the independence 
of total power on the flow rate in spot mode. From this 
comparison, the assumptions made in the cathode 
model appear to be better suited to the operation of 
larger and higher-current cathodes than originally 
intended. This suggests that the dominant physical 


processes for small hollow cathodes are different than 
large devices. 

Conclusions 

A model of the orificed hollow cathode, based on 
particle and energy balances, was presented. The 
model included the plasma in both the insert and orifice 
regions. In the orifice region, a single component, 
Maxwellian electron population was assumed to 
generate the ions and carry the bulk of the current. 
Models for the boundaries of the orifice plasma were 
presented. In the insert region, both a primary beam, 
thermionic and background Maxwellian electron 
populations were assumed to sustain the plasma and the 
current. The insert region current emission was 
constrained by the calculated current required by the 
orifice region. 

The sensitivity of the model to the various free 
parameters was assessed for low-current hollow 
cathodes. The cathode power consumption was found 
to be a strong function of the average excitation energy 
within the orifice, the fraction of excitation collisions 
resulting in ionization in the orifice, and the plasma 
potential in the insert region. Comparisons were made 
with experimental results to evaluate the utility of the 
model. A fundamental limitation of the energy balance 
calculation was discussed, and was found to be most 
critical for small-geometry devices at low-current. 
Additionally, the model when combined with 
experimental data suggested the presence of double- 
layers at the entrance and exit of the orifice. 
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Figure 1 - Illustration of the Orifice and Insert Model Approximations. 
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Figure 2 - Thermionic Current Density for Various Effective Work Functions. 
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Surface Temperature (K) 

Figure 3 - Dependence of the Thermionic Emission Current Density on the Plasma Potential Based on Double 

Sheath Analysis. Material Work Function = 2.0-eV. 
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Figure 4 - Electron Impact Cross-Sections for Xenon. 17 ' 18 ' 10 
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Figure 5 - Relative Dependency of Orifice Plasma Properties on the Fraction of Excitation Collisions 
Contributing to Ionization. AR6 Geometry at 0.50-A and 0.9-sccm. 
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Figure 6 - Relative Dependency of the Orifice Region Plasma Properties on the Average Excitation Energy. 

AR6 Geometry at 0.5-A and 0.9-sccm. 



f x in Insert Region (percent) 


Figure 7 - Dependency of the Insert Region Plasma Properties on the Fraction of Excitation Collisions 
Contributing to Ionization. AR6 Geometry at 0.50-A and 0.9-sccm. 



Figure 8 - Dependency of the Insert Region Plasma Properties on the Average Excitation Energy. AR6 

Geometry at 0.50-A and 0.9-sccm. 


NAS A/TM— 2002-2 11916 


13 






6 7 8 9 10 11 12 13 


Insert Region Plasma Potential (V) 

Figure 9 - Dependency of the Insert Region Plasma Parameters on the Plasma Potential. AR6 Geometry at 

0.50-A and 0.9-sccm. 



Insert Material Work Function (eV) 

Figure 10 - The Effect of the Material Work Function on the Predicted Insert Temperature for AR6 at 

0.50-A, and 0.9 seem. 


NAS A/TM— 2002-2 11916 


14 




•- 



u 


<u 

OX) 




0 ; 
OX) 


o 

•- 

o 

Ph 



Figure 11 - Sensitivity of the Predicted Power to the Free Parameters for a Hollow Cathode Operating at 

0.50-A. 



Figure 12 - Variation of Predicted Orifice Power Consumption Compared with Measured Cathode 
Performance. 22 The Calculated Dependency and Magnitude Indicate that the Bulk Averaged Plasma 

Properties Over-Predict the Ion Losses. 
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Figure 13 - Measured and Calculated Plasma Properties of a Hollow Cathode. 22 Calculated Volume 
Averaged Density in the Orifice Region Suggest the Presence of One or More Double Layers at the Upstream 

and/or Downstream Boundaries. 



Xe Flow Rate (seem) 

Figure 14 - Calculated Power Consumption in the Orifice of the Space Station Plasma Contactor Showing a 

More Modest Dependence Upon Flow Rate. 24 
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